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' Abstract 

We study central Z-boson production accompanied by rapidity gaps on either 
side as a way to gauge Higgs weak boson fusion production at the LHC. We analyse 
the possible backgrounds for the 6B-decay mode and show that these can be sub- 
stantially reduced. Special attention is paid to the evaluation of the gap survival 
factor, which is the major source of theoretical uncertainty in the rate of H, Z and 
W central production events with rapidity gaps. 



1 Introduction 



Hunting the Higgs boson(s) is the highest priority of the international high-energy physics 
programme. The Standard Model-hke Higgs boson should have a mass between the LEP2 
limit of 114 GeV and the upper bound of about 200 GeV, which is favoured by electroweak 
data [|l|. Within the MSSM, the light scalar Higgs boson is expected to be lighter than 
about 135 GeV, see for example 0. The focus now is on searching for the Higgs at present 
and forthcoming hadron colliders, namely the Tevatron and the LHC. 

To ascertain whether a Higgs signal can be seen, it is crucial to show that the back- 
ground does not overwhelm the signal. For instance, the major difficulty in observing 
inclusive production of the Higgs in the preferred mass range around 115 GeV via the 
dominant H ^ bb mode is the huge bb QCD background. In order to rescue the bb Higgs 
signal different options have been proposed in the literature. An attractive possibility to 
reduce the background is to study the central production of the Higgs in events with a 
large rapidity gap on either side, see for example [3-14]. An obvious advantage of the 
rapidity gap approach is the clean experimental signature - hadron free zones between the 
remnants of the incoming protons and the Higgs decay products. The cleanest situation 
is in the double-diffractive exclusive process: 

p^p ^p + H + '^p (1.1) 

where the plus sign denotes a large rapidity gap. However the cross section is expected to 
be rather small ^ , and as a consequence, the corresponding event rate appears to be 
too low at the Tevatron. Only at the LHC is there a chance of observing this exclusive 
Higgs production process [|TT], |T^, |TB|, [T^. Various effects cause a drastic reduction of 



the cross section for process ( [L . 1|) , for details see 0, |1^. First, the proton form factors 
strongly limit the available phase space in the transverse momentum of the produced 
Higgs, qt ~ ^/Rpi where Rp is the proton radius. Secondly, we have to account for the 
probability S'^ that the gaps survive the soft rescattering effects of spectator partons which 
may populate the gaps with secondary particles, see for example [|], 0, |2y]. Thirdly, the 
cross section is also suppressed by QCD Sudakov-like radiative effects 

The cross section is larger in the semi-inclusive case when the protons may dissociate, 

pp ^ X + H + Y (1.2) 

but the Higgs is still isolated by rapidity gaps. In this case there is no proton form factor 
suppression and the QCD "radiation damage" becomes weaker. Moreover, a significant 
contribution to process ( |1.2D comes from Higgs production via WW/ZZ fusion, i.e. qq — > 
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qqH. Since this process is mediated by colourless t-channel W/Z exchanges there is no 
corresponding gluon bremsstrahlung in the central region H, g I, 0, and thus Sudakov 
suppression of the rapidity gaps does not occur. Another characteristic feature of the 
vector boson fusion Higgs production process is that it is accompanied by energetic quark 
jets in the forward and backward directions. Recently, interest in this type of Higgs 



production process at the LHC has risen rapidly, see for example p!8|, p!9|, |2^, p3|. The 



particular importance of the electroweak fusion process is that it allows a determination 
of the Higgs coupling to vector bosons. It is worthwhile to note that the WW/ZZ fusion 
mechanism can provide a potential way to identify H bb decays at the LHC, if particular 
kinematic configurations with large rapidity gaps are selected, see for example ||18|| . 

The most delicate issue in calculating the cross section for processes with rapidity 
gaps concerns the soft survival factor^ S*^. This factor has been calculated in a number 



of models for various rapidity gap processes, see for example [19, 20, U4, Although 



there is reasonable agreement between these model expectations, it is always difficult to 
guarantee the precision of predictions which rely on soft physics. 

Fortunately, the calculations of S"^ can be checked experimentally by computing the 
event rate for a suitable calibrating reaction and comparing with the observed rate. As 
shown in 0, |21| the appropriate monitoring process for the double-diffractive mechanism 
is central dijet production with a rapidity gap on either side. To date, such a check 
has been the prediction of diffractive dijet production at the Tevatron in terms of the 
diffractive structure functions measured at HERA [^. The evaluation of the survival 
factor 5*^ based on the formalism of |1T9| , ^ appears to be in remarkable agreement with 
the CDF data (see also |p7| , p8|]). We expect that future measurements in run II of the 
Tevatron will provide us with further detailed information on S*^. 

As was pointed out in |2^, ^ , the survival factor for the gaps surrounding WW — > H 
fusion can be monitored experimentally by observing the closely related central produc- 
tion of a Z/W boson with the same rapidity gap and jet configuration. The discussion 
in 



29| , [30|] concerns the studies of Z/W + 2 forward jet production with subsequent lep- 



tonic Z decay in association with a rapidity gap trigger, which could allow a substantial 
suppression of the QCD-induced backgrounds. It is worthwhile to keep in mind that in 
different papers different criteria are used for the definition of the rapidity gaps. For ex- 
ample, in the approach of I^S], |30| no jets with px > 10 — 20 GeV are permitted within the 
gaps, while in [Rl p|, |11| the gap is required to be completely devoid of any soft hadrons. 



""^Recall that this factor is not universal, but is very sensitive to the spatial distribution of partons 
inside the colliding protons, which in turn results in the dependence on the particular hard subprocess 



as well as on the kinematical configurations of the parent reaction 19, pO| 
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Note that the determination of the gap survival factor in the vector boson mediated 
process is interesting in its own right, since here we can separate the contribution of the 
short transverse size component of the proton which so far has not attracted much 
attention, either theoretically or experimentally. 

The reader should be warned about the potential problems with the identification 
of rapidity gaps in the real life experimental environment at the LHC When the LHC 
operates at medium and high luminosity, the recorded events will be plagued by overlap 
interactions in the same bunch crossing (pile- up). However, as discussed in [|^, at least 
at the medium luminosity of 10^^ cm~^s~^ the gaps should be detectable. Using vertex 
reconstruction information, one can separate particles originating from the same vertex 
as the high-ET jets from those relatively low-p^ particles which arise from other vertices 
corresponding to the pile-up interactions. However, it is quite unlikely that this technique 
can be used for the super-LHC luminosity of 10^^ cm^^s^^. 

In this paper we develop the ideas of ||2^, ^ further by considering the decays of both 
(light) Higgs and Z bosons into bb pairs, the dominant decay channel of the former^. In 
each case we require two forward energetic jets, and rapidity gaps on either side of the 
centrally produced decay products. Both H and Z can be produced by electroweak vector 
boson fusion, for which gaps are 'natural', but the Z can also be produced via ^(a|) QCD 
processes, with both quarks and gluons exchanged in the t-channel, but as we shall see 
these have a smaller soft survival factor. Finally, there is a potentially large continuum 
bb background, which is again heavily suppressed when rapidity gaps are required. 

The paper is organised as follows. In Section ^ we discuss the calculation of the various 
signal and background processes at the parton level. In Section |] we perform detailed 
numerical calculations with realistic experimental cuts to determine the corresponding 
cross sections. In Section ^ we explore the consequences of hadronisation and estimate 
the gap survival probabilities, S"^, for the various processes under consideration. Finally, 
in Section ^ we combine the parton-level cross sections with the gap survival probabilities 
to give our final predictions for the cross sections. Section 1^ summarises our conclusions. 



2 Parton Level Calculation of Higgs, Z and W pro- 
duction 

In this section we assemble the various parton-level cross sections that are used in our 
analysis. We are particularly interested in the overall event rates for the various signal 
^For completeness, we also consider the + 2 jet production process. 
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and background processes, and in the kinematic distributions of the final-state parti- 
cles. All matrix elements used in the cross-section calculations are obtained using MAD- 
GRAPH [^. In all cases we work in the zero width approximation for the centrally- 
produced bosons. 

We begin by considering the fundamental signal process, ^(a^) Higgs production 
by WW, ZZ fusion: qq qqH (Fig. |I]). We assume that the Higgs is light, so that 
the dominant decay is into the bb final state. Because the momentum transfer is much 
smaller than the energy of the struck quark jets((pr) ~ Mi^/z), the jets are produced 
predominantly at small angle (i.e. large rapidity). Note that there is no exchange of 
colour in the t-channel, which leads to a suppression of hadronic radiation in the central 
region between the forward jets P, ||, 0. 

Representative Feynman diagrams for the analogous ^(a^) Z production process, 
qq qqZ and qq qqZ are shown in Fig. ^. They were first analysed in p9|, . Note 
that in addition to the WW fusion diagram. Fig. ^(a), the Z can also be radiated off either 
of the incoming or outgoing quark lines. Figs. 0(b) and (c). The characteristic topology 
of (b) is of a Z preferentially produced in the forward or backward region close in rapidity 
to one of the final-state quark jets. Requiring centrally produced Z decay products tends 
to suppress this contribution. Process (c) corresponds to s-channel production of the 
final-state qq pair, with the Z boson emitted off the incoming quark lines. It does not 
correspond to t-channel colour singlet exchange and is heavily kinematically suppressed 
by requiring a large rapidity separation between the jets. 

Similar remarks apply to W production. Representative Feynman diagrams for qq — > 
Wqq are shown in Fig. |^. Note that the central IV-production via 7 exchange correspond- 
ing to Fig. 



was recently discussed in |p2 




Figure 1: Higgs production via electroweak vector boson fusion. 



The above ^(a^) H and Z production processes both therefore give rise to rapidity 
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q,q q,q q,q q,q q Z 

(a) (b) (c) 

Figure 2: The three topologies for Zqq production via electroweak vector boson exchanges. 



q q q w q w 




Figure 3: The three topologies for W production. 



gap signatures between the forward jets and the central H and Z decay products. However 
there is a potentially important QCD ff^a^aw) background contribution to Z + 2 jet 
production where the internal electroweak gauge boson is replaced by a gluon. More 
generally, at this order indistinguishable background contributions can arise from any 
2 — >■ 2 scattering process (other than gg —>■ gg) where the Z is radiated off a quark 
line. Representative Feynman diagrams are shown in Fig. ^. By selecting forward jets 
and central Z bosons, in order to mimic the dominant Higgs configuration, the t-channel 
momentum transfer is minimised, and these QCD processes split into two types: t-channel 
quark (Figs. ^(a,b,c)) and gluon exchange (Figs. §(d)). Requiring rapidity gaps therefore 
suppresses both type of contribution, as will be discussed in Section ^ below. 

Finally, given that we are interested in the bb decay modes of both the Higgs and Z 
bosons, with two additional jets in the final state, there is a class of ^(a|) pure-QCD 
background processes of the form ab cd + bb with a...d = q,g, examples of which are 
shown in Fig. ^. We will consider the corresponding cross sections in the following section, 
with the additional requirement that mj,^ ~ Mz- 
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3 Signal and Background Rates and Properties 

In this section we calculate cross sections and kinematic distributions for the processes 
described in the previous section, for proton-proton collisions at the LHC using a repre- 
sentative set of cuts on the final state particles. We have in mind final states with a jet 
registered in a forward detector with pseudorapidity rji > r^min, another produced back- 
wards with ri2 < — ?7min, and the H, Z and W decay products produced centrally, with 
rapidity \yH,z,w\ < Umax- We impose a minimum transverse momentum cut pTmin on all 
final-state jets. 



3.1 Total cross sections 

Figure ^ shows the total cross section for Higgs, electroweak Z and W, and QCD Z 
production (with no branching ratios included) as a function of a cut on the minimum 
transverse jet momentum prmin- The Higgs mass is Mh = 115 GeV and the leading-order 



MRST98LO parton distribution set is used. Note that only for H production is the 
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cross section finite in the limit PTmin ^ 00. In addition, the possibihty that the final state 
jets in Z and W production originate in the splitting process 5'*,7* qq (for example, 
see Fig. 0(c)) requires a jet separation cut. The minimal way to do this is simply to 
require that one of the jets is produced in the forward hemisphere and the other in the 
backward hemisphere, i.e. rji ■ ri2 < 0. When we come to consider 'realistic' cuts, in 
particular to isolate the jets from each other and the H and Z decay products, we will 
impose a large rapidity separation cut in which one jet is produced far forward and one 
far backward: |?7i|, \ri2\ > rjmin, Vi ' V2 < 0. For the Higgs production process, which has 
no infrared or coUinear singularities, the imposition of prmm and rji ■ ri2 < acceptance 
cuts simply reduces the cross section slightly (by approximately 25% for a broad range of 
PTmin values), see Fig. |. 

Figure ^ shows that there is a strong ordering of the cross sections cr(Z, QCD) ^ 
a{Z,EW) ^ (^{H), with a{Z,QCD) exhibiting the strongest dependence on pTmin- The 
W cross section has a stronger infra-red singularity as pTmin than the corresponding 
Z cross section, due to the soft photon singularity present in the extra diagram with 
respect to the Z production process involving the triple gauge boson vertex (Fig. 0(a)). 
This is shown more clearly in W/Z cross section ratio plot. Fig. ^. The Higgs cross section 
is only weakly dependent on the mass Mh, decreasing by a factor of 2 as Mh increases 
from 100 GeV to 200 GeV, see Fig. |. 

Note that all the above cross sections are evaluated in the zero Z/W width approxi- 
mation and at leading order in perturbation theory. In particular, in the QCD Z + 2 jet 
calculation the scale of the strong coupling is not determined, and there is a non- 
negligible scale dependence uncertainty as a result. We use as = as^M^) throughout. 

■^The possibility of exchanging a massless photon or gluon in the t channel gives rise to an infrared 
singularity in the electroweak and QCD Z + 2 jet production processes as pTmin ~^ 0, see Figs. || and ^. 
For exclusive pp ^ p + X + p collisions this singularity is cut off by the i,nin-effect. 
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Figure 6: Total cross sections for (H, Z, W) +2 jet production in pp collisions at ^/s = 14 
TeV as a function of the cut on the jet transverse momentum, pTmin- Rapidity cuts 
on the final state jets are also imposed, as indicated. 



One could also argue for a smaller scale characteristic of the transverse momenta of the 
forward jets, e.g. as = CisiPTmm)- ^^^^ discuss the impact of such a choice on our 
predicted event rates in Section H. 



3.2 Distributions 



Our objective is to find a set of selection cuts that minimises the background while not 
affecting drastically the Higgs, Z and W rapidity gap signal. We begin by calculating the 
transverse momentum and rapidity distributions of the jets in qqH production, Figs. || 
and 0. Evidently the jets are predominantly produced with transverse momenta of order 
M\y/2 ~ 40 GeV, and with a rapidity separation of around 5, see Fig. |TT|. Notice the 
small excess around Ar] ~ 1/2. This is caused by the contributing process qq Hqq in 
which rrijj ~ Mz, i.e. the Higgs is produced in association with a Z (or W) boson which 
subsequently decays into a qq pair, see Fig. |12[ This is more clearly seen in the dijet mass 
distribution. Fig. |l^. 

Requiring the jets to be well-separated in rapidity forces rrijj to be large and this 



resonant contribution is strongly suppressed. For example. Fig. 13 also shows the dijet 
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Figure 7: The ratio between tiie total cross sections of the W and Z (electroweak) production 
processes as defined in the previous figure. The W cross section is 'more divergent' 
than that for Z production at low pT, because of the extra photon-exchange diagram 
involving the triple gauge boson vertex (Fig. ^(a)). 
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Figure 8: The ^/s = 14 TeV qq qqH cross section as a function of the Higgs mass Mh- 
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Figure 9: Jet transverse momentum distribution for qq qqH. The peak is at around 
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Figure 10: Rapidity distribution of jets for qq qqH with pTmin = 40GeV and r]i.rj2 < 0. 
The vertical scale is normalised arbitrarily. 
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Figure 11: Dijet rapidity difference for qq — > qqH (with 7/1.772 < 0) as a function of the PTmin 
transverse momentum cut. The vertical scale is normalised arbitrarily. The gap 
narrows marginally as the cut is raised, as expected from kinematics. The small 
excess seen at low Ar/ is discussed in the text. 



mass distribution for |?7i 2I > "/^min = 2. 

The jet rapidity distribution for electroweak qqZ production is shown in Fig. 



Comparing with Fig. [ly for qqH , we see that the jets produced with a Z are more 
uniform in rapidity, The Wl^-fusion' diagrams of Fig. ^(a) stih produce jets with a large 
separation, but the central region is now filled in by contributions from the other non- 
fusion 'Z-bremsstrahlung' processes. Figs. @(b,c). Electroweak W production has very 
similar characteristics to electroweak Z production. 

For the QCD background to electroweak Z production, the jets are produced much 
more centrally, see Fig. Requiring a jet in each forward/backward hemisphere leads 
to a typical rapidity separation of about 3, as shown in Fig. |16|, which is significantly less 
than for either H or electroweak Z production. There is no natural rapidity gap, as for 
the t-channel colour-singlet exchange processes. 



3.3 Selection cuts 

We can now proceed to define a set of selection cuts that leads to a sample of H, Z and W 
events with the potential to exhibit rapidity gaps. Since our primary goal is to calibrate 
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Figure 12: Contribution to the ^{a'^r) electroweak process qq — > Hqq that resonates when 
rrijj ~ Mz- 



the gap survival for Higgs production, we will concentrate first on the bb decays of H and 
Z, the latter produced either via electroweak or QCD processes. 

When considering the bb decay modes of both the Higgs and Z bosons, we must include 
also the important irreducible background from QCD ^(al) 66 + 2 jet production, see 
Fig. ||. Such processes give a continuous distribution of m^j masses, and in what follows 
we impose a cut of jm^^ — Mz\ < 10 GeV to select those background events that mimic 
Z ^ bb decay. 

The configuration we have in mind has one jet registered in a forward detector with 
V > Vmin, another produced backwards with r] < —rj^in, and the two 6 jets from H and Z 
decay produced centrally. From the results of the previous section, such a selection will 
in principle preserve the bulk of the Higgs signal while suppressing the (non-gap) QCD 
Z and 66 + 2 jet production. 

For both ATLAS and CMS, the forward hadron calorimeters cover approximately 
3 < Ir^l < 5, and so we will require our forward dijets to be produced in this region of 
rapidity, with pt > PTmin = 40 GeV. In order to separate the H, Z decay jets from the 
forward jets, we require |?7f,| < 1.5, and prt > 10 GeV0. Although these cuts are designed 
to reflect the 'natural' characteristics of qqH production, they do result in a non-negligible 
loss of signal rate, even before 6-tagging efficiencies etc. are taken into account. This is 
illustrated in Table |I|, which quantifies the effect on the cross section of applying the cuts 
sequentially. One can see that imposing forward jet cuts has the largest impact on the 
cross section, and indeed this is the case for all the processes considered. 

Figures ^ and show the cross sections at ^/s = 14 TeV as a function of pTmin for 

^The typical transverse momentum of the jets in both tlie signal and background processes is ~ Mz/2, 
and this cut does not have any significant effect on the event rates, see for example Table |^ below. 
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Figure 13: The dijet invariant mass for qq — > qqH shows a double resonance around mz and 



mw ~ due to diagrams such as Fig. 12. The vertical scale is normalised arbitrarily. 



When the jets are required to be separated and forward in rapidity the effect is 
irrelevant. 
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Figure 14: Rapidity distribution of jets for electroweak Z production. The vertical scale is 



normalised arbitrarily. The gap is narrower than for the Higgs signal (Fig. IC). A 
jet transverse momentum cut of prmin = 40 GeV is applied, as is rji.r]2 < 0. 
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Figure 15: Rapidity distribution of jets for QCD Z + 2 jet production. The vertical scale 
is normalised arbitrarily. A jet transverse momentum cut of PTmin = 40 GeV is 
applied, as is 771.772 < 0. 




Figure 16: Comparison of the dijet rapidity differences for the electroweak and QCD Z pro- 
duction processes. The vertical scale is normalised arbitrarily. A jet transverse 
momentum cut of prmin = 40 GeV is applied, as is ?7i.?72 < 0. 
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Cut Imposed 


Cross Section for qq — > qqH at px > 40 GeV 


% of Initial Cross Section 






inn% 


VirfH hh) 

J— *1 \J-J- ~ t/t/ 1 




71 9% 


T^^ • T7o <' n 
'/I '/2 ^ 


2.47pb 


50 8% 


Ar^j > 6 


0.495pb 


10.2% 


l^.l >3 


0.0990pb 


2.04% 


< 1.5 


0.0465pb 


0.957% 


VTh > 10 GeV 


0.0463pb 


0.953% 



Table 1: Loss of qq qqH cross section at y/s = 14 TeV with Mh = 115 GeV in applying 
selection cuts and the bb branching ratio. 



all processes. The Higgs production cross section is reduced by a factor of ~ 100 and 
the electroweak Z production by ~ 1000 in comparison with Figure ^ The cuts reduce 
the Z production QCD background by a factor of ~ 10000. As already mentioned, in 
evaluating the pure QCD bb production cross sections we further impose the restriction 
that the dijet invariant mass be within 10 GeV of Mz- 



4 Gap Survival Probability 
4.1 Parton Level 

In the previous sections we have presented cross sections for Higgs, Z and QCD bb pro- 
duction processes for events with rapidity gaps at the parton level. We take the definition 
of rapidity gap to mean that there should be no minijets with a large {pt > 10 GeV) 
transverse momentum within the gap region. As discussed in the introduction, the selec- 
tion of rapidity gap events improves the signal to background ratio because gaps are a 
characteristic feature of the vector boson fusion process, whereas they are not for QCD 
Z and continuum bb production. Of course, the results presented in the previous section 
should be corrected to account for the rescattering of spectator partons, that is the pos- 
sibility that another pair of initial, fast partons interacts inelastically in the same event. 
Secondaries produced in this inelastic interaction may fill the gap and the probability, 
S"^, for the gap to survive depends on the criteria used to select the gap. Insofar as we 



^This is because of the difterence in rapidity distributions in the H (Figs. 10 and |ri|) and Z (Figs, 
and |l^) cases which is caused by the process shown in Figs ||(b),(c) and Fig. ^(d) where the quark jets 
are closer to each other. 
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Cuts: Tij.ri^ < ; It], ,1 > 3 ; ItiJ < 1.5 ; p.^^ > 10 GeV 
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Figure 17: Cross sections at ^/s = 14 TeV for Higgs production processes after application of 
the cuts described in the text. 



Cuts: rij.ri^ < ; Irj, ,1 > 3 ; IrjJ < 1.5 ; p.^^ > 10 GeV 



'1,2 



Tb ■ 



le+04 
1000 
100 
10 

1 

0.1 
0.01 
0.001 
0.0001 







- \ 

\ 








qq qqZ (Electroweak); Z ^ bb 

- ~ - qq, qg, gg ^ Z + 2 jet (QCD); Z ^ bb 

qq ^ qqbb; m^^ - ± 10 GeV 

ffp p?bb: m. _ = m_ + 10 GeV 


\ 

= \ ^ 

- N 

* s 


V _ 

' 






-\ 
- \ 




■ — ■ > 






c 


ig ^qg 


bb 


z 

= m^± 


10Ge\ 


7 


r-\ 














































* — . 


































1 








1 








1 





10 20 



30 



40 
Pt 

'^T mm 



50 60 

(GeV) 



70 80 90 100 



Figure 18: Cross sections at ^/s = 14 TeV for Z production processes after apphcation of the 
cuts described in the text. 
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require only to have no high pt particles (or minijets) within the gap interval, the effect 
is not too large. The probability to produce high secondaries is relatively small and 
the corresponding survival factor 5*^ ~ 0.7 — 0.9 (depending on the pt cut) is more or less 
close to one, see for example p3[. For any specific kinematics (and pt cuts), the value of 
5*^ for such 'parton level' gaps may be estimated using, for example, one of the options of 



the PYTHIA Monte Carlo program |]3^, or an 'exponentiation' model such as that used 



m 



4.2 Hadron Level 



As seen in Figs. |T^, |18|, the QCD-induced bb background is still large. It exceeds by two 
orders of magnitude the Z/H cross sections and it is therefore necessary to suppress the 
background further. This can be done by requiring a completely clean gap, i.e. without 
any soft hadrons. Indeed, all the QCD processes we consider are characterised by gluon 
(or quark) t-channel exchange, which unavoidably produces a colour flow along the gap. 
During hadronisation this colour flow, in turn, creates plenty of soft secondaries which fill 
the gap. On the other hand, there is no such effect for the electroweak graphs (Figs. |1], 
l^a, Pd) since the vector boson exchange is colourless. This means that if we require clean 
hadron-level gaps we can immediately discard the diagrams of Fig. |^. The only way to 
create a gap in a QCD induced event is to screen the colour flow (across the gap) by an 
additional gluon (or quark) exchange; that is, to consider graphs of the type shown in 



Fig. 19 



Pti 




PT2 




Pti 



PT2 



Figure 19: Screening of QCD dijet + bb production via gluon exchange. 

Note that in leading order we can screen the colour flow in both gaps (above and 
below the 66-pair) with only one additional t-channel gluon, with momentum Qt say. 

^ Note that the large survival probability ~ 0.8 used in ||2^] corresponds just to parton level gaps, 
and was calculated using the model of . 
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The price we pay for this screening is a factor of as supplemented by the dQ'^ loop 
integration in each amplitude; that is {j as ■■ ■ d'^QrY ^^e cross section. At first sight, 
the major contribution comes from the small Qt region where the QCD coupling is 
larger. Moreover, the integral takes the form 



(4.1) 



and has infrared logarithmic divergence at Qt <^ Prjeti,2- However, this divergence is cut 
off by the effects of higher order double logarithmic QCD radiation, see for example p, ^ 



21[| . The point is that a small Qt gluon screens the hard gluon at rather large distances 
r ~ 1/Qt only. Thus a 'hard' gluon gi=i,2 may emit a new 'semihard' gluon jet, with 
transverse energy Et ranging from Qt up to g^T = Wxi — Qt\ in the whole rapidity gap 
interval Ar^j. The leading logarithms come from the Qt <^ Et <^ PTjet domain where the 
expected mean number of these secondary gluons is 

fi^^^A,,.,,,^. (4.2) 

TT Q^ 

At the amplitude level the corresponding suppression factor describing the probability 
for not having such an emission (which otherwise destroys the gap) has the Sudakov-like 
form 

NrCt Q A 



exp(-n,/2) = ( ^ ) . (4.3) 



Including this factor in the loop integral, we eliminate the infrared divergence and obtain 
the probability. Pa (a = q(l,qg,gg depending on the initial state), to screen out the octet 
(gluon-like) colour flow in qq {gg or qg) interactions. 



Pa = Ca{ / «5(QT)^exp ^ / as{Q")^ = a 



Q't { 27r Jq^ ' Q'^ }J \NAV 



n, 

(4.4) 

Here A?7 = Ar/i + Ar72 is the overall length of the gaps and, within leading logarithm 
accuracy, we have put the upper limits in the Qt {Q') integration equal to the minimum 
Pt of the jets. In order to arrive at the right-hand side of Eq. ( [4. 41) it is convenient to 
recast the integral in Eq. ( [4.4| ) as 



27r 



iVcAr? 



dJe^-^{-J{pT^in,QT)) (4.5) 



with 
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Performing the integration we neglect the term exp(— j7'(pTmin, Qo)) corresponding to the 
lower limit of integration. This can always be done safely if we can continue the pertur- 
bative calculation down to the (rather low) scale where the quantity ^^(Qo) ' becomes 
large. Instead of the conventional double logarithm expressions (Eqs. ( [4.2|) and ( [4.3|) ) 
with a fixed coupling as, in Eq. (|4.4|) we have used the running coupling in order to 
demonstrate that the result does not depend on whether one accounts for the running as 
or not. The colour factors Ca are 

C| _Cf _2 _ N! _9 _ CpN, _ 1 

(iV2-l) 2iV, 9' 8' 2' ^ ' 

A more precise way to calculate the contributions of Fig. |19| including QCD radia- 
tive effects is to replace the two gluon t-channel exchange by the non-forward BFKL 
amplitude |]36[. For the asymmetric {Qt ^ qu) configuration the non-forward ampli- 
tude contains the double logarithmic factor of Eq. ([4.3|) , while the single logarithmic 
(~ 0{asAri)) contribution in this asymmetric kinematical situation is suppressed, giving 
a less than 10% correction to the amplitude (see ^1, ^ for a more detailed discussion). 



Thus we come back to the result of Eq. ( |4.4| ). Strictly speaking, besides the suppression 
factor Eq. (|4.3|) hidden in the BFKL amplitude, there should be another Sudakov-like 
double logarithmic form factor which reflects the absence of QCD radiation in the inter- 
val of gluon transverse momentum between prjet and half of the boson (or bb) mass, M/2. 
However, in our case the transverse momentum of the jets is px > Prmm = 40 GeV, which 
is close to half the boson mass Mz^h/'^- Therefore the form factor becomes close to one 
and we can neglect it. 

Another point we have to take into account is the fact that now the fe6-pair may be 
produced in a colour singlet state only, and the ordinary gg bb hard subprocess cross 
section (which includes both colour singlet and octet contributions) should be replaced 
by the pure colour singlet cross section p8 



rfa''^^', PP . nal 1 \( ^ 1El\ ( _ 2m2\ 



-{gg 



N^-l dt ''''' {N^ -1)E^M^6 



(4.; 



where (3 = y 1 — -j^ and is the quark mass. Note that for the colour singlet production 
case there is an additional colour factor l/{Nl — l) which suppresses the QCD background, 
as the two colliding gluons are forced to have the same colour. 



4.3 Quark Exchange 

It is more difficult to screen the colour triplet flow originated by the quark exchange 
which we deal with in the electroweak and QCD Z + 2 jet processes shown in Figs. 0c, 
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|c or ^a,b,c. For example, to screen the quark colour in Fig. | we have to replace the 



graphs Figs. §a,b,c by those of Fig. |2^. Due to the spin 1/2 nature of a quark, the large 





m 




(a) 



(b) 



(c) 



Figure 20: Screening of QCD dijet + Z production via quark exchange. 



rapidity gaps are suppressed at the amplitude level (in comparison with the corresponding 
Figs. ^,b,c amplitude contribution) by the factor e~^'^^'^ (i.e. a factor e"^'' in the cross 
section). On the other hand, it is known that the loop with two t-channel fermions may 



contain a double logarithm (see p9| , 0). One logarithm comes from the transverse (Qt) 
integration, while another logarithm (in the real part of the amplitude) originates from 
the dm? /rn? integral over the (virtual) mass of the upper s-channel particle in the loop 
(assuming that the contour of the Feynman integral is closed on the pole corresponding 
to the lower s-channel particle). In our kinematics, where a Z boson is emitted in the 
centre of the rapidity gap interval, we obtain a logarithm when the mass, m, runs from 
m? = maxjQl^, a/sMI^} up to m? = s (here s is the incoming parton energy squared and 
Ml J, = Ml + \ptz\)- That is, the mass integral gives / ^ < The Qt integration 
does not give a logarithm in the case of Fig. pO|a, but for the amplitudes corresponding to 
Fig. [20|b,c, a logarithmic integral appears in the domain pf-jct ^ Q% ^ Thus from 

the Fig. |20|b,c loop integrals we may expect a ^{Aiff enhancement. However, with our 
large rapidity gap (Ar/ ~ 6) the whole factor - [^li^rjfe-^'^'^f = 0.45(^)2 ~ lO"^ is 
very small. Besides this, after the parton-level cuts described in Section p.2| are applied, 
the original parton-level contribution of the diagrams with a (t-channel) quark exchange 
is strongly suppressed. Therefore we neglect these contributions. 



4.4 Soft Survival Probability 

Returning to our original processes, we keep now only the graphs with either vector boson 
or two gluon (Fig. |T9|) exchange across the rapidity gaps, and multiply the corresponding 
parton-level cross sections by the appropriate gap survival probability 5*^. However, as 
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we now require there to be no hadrons (even with a rather low px) in the gap interval, 
we have to take account of any soft interactions of the spectator partons. 

Instead of using a Monte Carlo simulation, it is better to choose a model based on 
Regge (Pomeron) theory tuned to describe soft interaction data at high energies. We 



will use the model of Ref. This is based on the two-channel eikonal formalism, 

which reproduces all the main features of the soft {(Jtot, do'ei/dt) cross section data in the 
ISR-Tevatron energy range. Recall that the two channels of the eikonal correspond to 
two eigenstates which have different absorptive cross sections (i.e. different rescattering 
probabilities). Assuming the same (momentum and spatial) distributions of quarks and 
gluons in both components of the incoming proton wavefunction (that is, in both eigen- 
states of the eikonal) the model predicts for all our processes S*^ = 0.1 at LHC energies. 
In other words, by requiring gaps at the hadron level we decrease the cross section by 
an order of magnitude. At first sight, the gap survival probability S"^ = 0.1 reflects the 
rescattering of soft spectator partons and should, therefore, be universal for any process 
which has a gapQ. However, this is not completely true. First, the value of depends 
on the spatial distribution of parton spectators and therefore on the characteristic impact 
parameter (6^) difference between the two colliding protons |jl9|, For example, in the 
case of exclusive Higgs boson production, pp p + H + p via photon-photon fusion, the 
transverse momenta of the photons are very small. Hence the impact parameter Bt is very 
large. The probability of soft rescattering in such a highly peripheral collision is small, 
and the value of (~ 0.9) is close to one |Tl|. Secondly, there is a difference in the 
momentum distributions of partons in a different (diffractive eigenstate) component of 
the incoming proton wavefunction; it is reasonable to expect that the component with a 
smaller cross section contains more valence quarks (and 'hard' large-x partons), whereas 
the component with a larger cross section has more low-x gluons. This possibility was dis- 
cussed in In such an approach, the model describes the breakdown of factorisation, 
in that there is about a factor 10 difference between the 'effective' Pomeron structure 
functions measured in diffractive deeply inelastic interactions at HERA and diffractive 
high-ET dijet hadroproduction at the Tevatron ||26[]P|. 

In the present context, as the background 66-pairs are produced predominantly in 
gluon-gluon collisions, the gap survival probability for the QCD background is a little 
smaller than for Z{H)-hoson production via vector boson fusion where we deal with 
incoming quark-quark interactions (see Figs. and 0). Using the formalism of Ref. ^0 



^The only difference may be caused by the Sudakov-like form factor tliat accounts for tlie absence of 
QCD gluon bremsstralilung in a specific liard subprocess. 

^The difference is explained simply by the fact that the gap survival factor is 5^ ~ 0.1 for proton- 
antiproton collisions, whereas 5^ ~ 1 in deep inelastic scattering. 
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we obtain for the kinematics (cuts) described in Section ^ 



Si = 0.31; Si = 0.31; Sl^^^j, = 0.27. (4.9) 



These survival factors are much larger than in the original model |T9[ because for the 
case considered here, of large rapidity gaps and large jet transverse momenta, we select 
mainly fast incoming partons and valence quarks which belong to the second component 
of the proton wavefunction. This component has a smaller absorptive cross section^ In 
this case the QCD background is additionally suppressed 2.5 times. Note that both 



versions of the model describe the diffractive dijet CDF data |^ well enough. On 
the other hand, in processes with large rapidity gaps at the LHC the uncertainty in the 
soft survival factor S*^ may be rather large. It will therefore be important to study such a 
process experimentally. A promising way to study the survival probability 5*^ in different 
components of the incoming proton wave function (i.e. the dependence of 5*^ on the PTjet 
and rapidity cuts) is to measure QCD dijet production with rapidity gaps on either side 
of the dijet pair. Here the cross section is much larger (especially for gluon-gluon induced 
dijets) and it is easy to study the gap survival factor S"^ under the various kinematic 
conditions: px of the fast (large rf) jets, size of the rapidity gaps, dijet mass, etc. In this 
way we can emphasise the role of the incoming valence quarks, sea quarks or gluons in 
different x and scale /i^ ~ p\ domains, and hence choose the configuration where one or 
other component of the wavefunction dominates. 

Note that, depending on the jet-finding algorithm, some soft hadrons may or may not 
be attributed to a particular fe-jet. Therefore, one has to be more specific in the definition 
of the rapidity gap on the hadronic level in the presence of the high-p^ jets. It looks 
plausible to select the gap by the requirement not to have hadrons within the gap range, 
apart from the cones of a fixed size Ai? ~ 1 around the jet directions. In a real life 
experiment, jet-finding algorithms should be utilised in optimising the value of Ai?. Soft 
survival factors S"^ are practically independent of the Ai? value at Ai? < 2. 



5 Results 

Fig. |21| shows the cross sections after hadronisation for central production of Higgs or 
Z with rapidity gaps and subsequent decay to bb as a function of PTmin of the proton 



^Under the extreme hypothesis that all valence quarks belong to the second (low Gabs) component 
while gluons and sea quarks are concentrated in the first component (with a larger cross section) we get 

Sl = Sl = 0.26 and 5^ ^^^^ = 0.10. 
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Cuts: rij.ri^ < ; Irj^ ^1 > 3 ; \\ < 1.5 ; p^^ > 10 GeV 



- ^ 1 
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pp ^ X + H + Y (nijj = 1 15 GeV) ; H ^ bb 

pp^X + Z + Y;Z^bb 
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. . . . pp ^ X + bb+ Y (QCD) ; m^-^= ± 10 GeV 
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Figure 21: Hadron-level cross sections at ^/s = 14 TeV for inclusive Higgs and Z production 
with subsequent decay to bb and their respective QCD bb backgrounds. Cuts are 
applied at the parton level as discussed in the text. 



remnant jets. It also shows the expected background of QCD bb events that display 
the same kinematical configuration. These are calculated using as a starting point the 
parton level cross sections after application of cuts, namely Figure |l^ for Higgs production 
and Figure |18| for Z production. The QCD-induced cross sections (both the QCD Z 
production of Fig. ^ and the direct QCD bb production of Fig. |^) are then multiplied 
by the probability to screen out the colour octet contribution for the relevant initial state 
of either qq, qg or gg (Eq. ([4.4|)). To take into account the fact that the bb pair in 
the background processes can only be produced in the colour singlet state the ordinary 
gg bb cross section is replaced by the pure singlet cross section, Eq. (|4.8|) . Finally both 
the signals and backgrounds are multiplied by the relevant soft survival probability of 
Eq. il). 

We see that, as long as we stay away from the low pTmin region, the signal for Higgs 
production is comparable with the QCD background, even exceeding it above PTmin = 
50 GeV. The cross section for Z production with rapidity gaps is less than that for Higgs 
production over most of the plot. This is because the branching fraction to bb is much 
lower than for the Higgs. Exceptions to this occur at low pTmm where the effect of the 
infrared singularity makes its presence felt and at very high pTmin, explained by the fact 
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that the parton-level Higgs cross section falls more rapidly, as shown in Fig. ^ The 
backgrounds show an extremely strong dependence on pTmin, falling by five orders of 
magnitude as one varies PTmm from 2 GeV to 100 GeV. The QCD bb background with 
the invariant mass of the bb pair taken around the Z mass is approximately 80% higher 
than that evaluated around the Higgs mass. 

It should be noted that in these calculations we have taken as = as{M^). An ar- 
gument could also be made that the characteristic scale should be that of the transverse 
momenta of the forward jets, i.e. as = ctsiPTmm)- This would affect the back- 
grounds in such a way as to increase them by approximately 30% if we take the typical 
Pt to be 40 GeV. 

Let us emphasise that up to now we have not addressed the experimental issues. In 



particular, the predictions given in Fig. |2T| should be modified to account for the 6-tagging 
efficiency Sb- This, in turn, is correlated with the probability V{g,q/b) to misidentify a 
gluon (or a light quark) jet as a 6-jet. Recall that the rate of the gg-dijets exceeds the 
66-yield by two orders of magnitude. As discussed in it is feasible to expect for the 
two 6-jets {ebf = 0.6 with V{g,q/b) = 0.01. 



6 Conclusions 

The weak boson fusion mechanism can provide a promising way to detect a light Higgs 
boson at the LHC, see for example ||18|, |2^. The selection of events with large rapidity 
gaps and energetic large px (quark) jets in the forward and backward directions allows 
the suppression of the QCD bb background down to a level comparable to the signal. 
Therefore, it becomes feasible to observe a light Higgs boson via its main H ^ bb decay 
mode in addition to the usually discussed rr and WW* channels, see for example |l^ . 



The cross section for the production of a 115 GeV Higgs boson in association with 
rapidity gaps at the LHC is expected to be about 15 fb (for pr > 40 GeV). Therefore, 
for an integrated luminosity of 30 fb~^ planned for the first two or three years of LHC 
running, about 400 events can be observed. 

Note that our cuts were not finally optimised for the particular ATLAS/CMS con- 
ditions. Thus, the significance of the signal may be improved by allowing asymmetric 
configurations with some minimal Arj between the high-p^- jets instead of the require- 
ment 1 771,2 1 > 3, i]i ■ ri2 < 0. Such a kinematical choice was considered, for instance, in 



Ref. |18|. It is shown that this will noticeably improve the significance of the signal. 



An important ingredient in the evaluation of both the signal and the background in 
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the bb + 2 forward jet events is the soft survival factor 5*^, defining the probabihty that 
the gaps survive the soft pp-scattering. Recall that though this factor can be computed 
within the framework of existing models for soft rescattering, it is always unwise to rely 
on the precision of models based on soft physics. Fortunately, the soft survival factor for 
the gaps surrounding WW — > H fusion can be monitored experimentally by observing 
the closely related central production of a Z boson with the same rapidity gap and jet 
configuration [p^, IST 



As was emphasised in ||30| , ^ , the rr and WW* decay channels with rapidity gap 
kinematics can give a rather high significance for the observation of a light Higgs. In the 
TT case the main background results from the tail in the rr mass distribution generated by 
the Z —>■ T^T^ decay. Again, the experimental observation of Z boson central production 
allows one to control and monitor such a background. 

It is worthwhile to mention that the experimental determination of the gap survival 
factor in the processes under consideration is interesting in its own right, since it pro- 
vides important information on the incoming proton wavefunction. Note that since the 
incoming partons in the subprocess qq ^ q + (bb) + q are rather hard, the factor S"^ 
depends on the model assumptions more sensitively than, for example, in the exclusive 
diffractive production case pp ^ p + bb + p, see JTT], . As was demonstrated in Sec- 
tion |4.4| (see footnote^) the results strongly depend on how the partons in the proton are 
distributed between the different diffractive eigenstates. Currently our information on 
these distributions is rather limited. 

This paper concentrates on the detailed analysis of central Z boson production ac- 
companied by rapidity gaps on either side and two forward jets at the LHC The QCD 
background processes for Z + 2 jet production in the rapidity gap environment are ad- 
dressed in detail. We evaluate the soft survival factors S'^ for various processes under 
consideration. Finally, we note that it will be important to extend our work by incor- 
porating a realistic Monte Carlo simulation, which will allow detector simulation to be 
included. We believe that the results presented in this paper make such an effort worth- 
while. 
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